Continuum approach to wide shear zones in quasi-static granular matter 
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Slow and dense granular flows often exhibit narrow shear bands, making them ill-suited for a 
continuum description. However, smooth granular flows have been shown to occur in specific ge- 
ometries such as linear shear in the absence of gravity, slow inclined plane flows and, recently, flows 
in split-bottom Couette geometries. The wide shear regions in these systems should be amenable to 
a continuum description, and the theoretical challenge lies in finding constitutive relations between 
the internal stresses and the flow field. We propose a set of testable constitutive assumptions, in- 
cluding rate-independence, and investigate the additional restrictions on the constitutive relations 
imposed by the flow geometries. The wide shear layers in the highly symmetric linear shear and 
inclined plane flows are consistent with the simple constitutive assumption that, in analogy with 
solid friction, the effective-friction coefficient (ratio between shear and normal stresses) is a constant. 
However, this standard picture of granular flows is shown to be inconsistent with flows in the less 
symmetric split-bottom geometry — here the effective friction coefficient must vary throughout the 
shear zone, or else the shear zone localizes. We suggest that a subtle dependence of the effective- 
friction coefficient on the orientation of the sliding layers with respect to the bulk force is crucial 
for the understanding of slow granular flows. 



T3 

O 
o 



> 

in 
o 

in 
o 

-I— » 



i 

-a 
c 

o 
o 



X 



I. INTRODUCTION 

Granular flows show a very wide variety of behaviors, 
and while the microscopic dynamics of dry cohesion-less 
grains is simple and well understood, there is no gen- 
eral theory describing their emergent macroscopic prop- 
erties. Flowing grains can roughly be classified into three 
regimes by the relative importance of inertial effects [1] . 
For strong external driving the grains form a gaseous 
state. Here particle interactions are dominated by binary 
collisions, and this regime is well captured by modified 
kinetic theories [2-5]. On lowering the driving strength 
the flow becomes denser, with collisions becoming cor- 
related and often involving several particles at once. In 
this regime inertia is still important, but kinetic theo- 
ries become increasingly difficult to justify and apply [6]. 
On further lowering the driving strength, the granular 
media enter the quasi-static regime where inertial effects 
are negligible. The grains form enduring contacts, lead- 
ing to highly complex contact and force networks. The 
modeling of these flows is still in its infancy, and there is 
no general approach which, for given geometry and grain 
properties, predicts the ensuing flow fields. 

Most slow granular flows can not be considered 
smooth. Their flow fields vary strongly on the grain scale. 
For example, in many experimental realizations one ob- 
serves shear localization where the flow of the material 
is concentrated into a very narrow shear band [7-10]. 
In such situations the flow can be modeled as two solid 
blocks sliding past each other. If the shear stress, t, is 
simply proportional to the normal pressure, P, then these 
materials are referred to as ideal cohesion-less Coulomb 
materials. Many formulations of granular flow focus im- 
mediately on this narrow shear band regime [8, 10]. 



However, there are a number of systems that dis- 
play smooth velocity fields and wide shear zones. These 
should be amenable to a continuum description. Among 
these is the planar-shear cell without gravity, which has 
been examined numerically in [11-16] (Fig. 1 a). An- 
other example is slow flows down inclined planes, which 
in simulations of 3D systems appear to reach a quasi- 
static state [17] (Fig. 1 b). Though conceptually simple, 
both of these situations are hard to realize experimen- 
tally. In recent experiments, using a modified Taylor- 
Couette cell with split bottom (Fig. 1 c), robust and 
wide shear zones were obtained in the quasi-static, dense 
regime [18-20]. While we do not expect there to be a 
" universal" continuum theory of granular flow, these ob- 
servations strongly suggests that there is a continuum 
theory with its own domain of validity, that should cap- 
ture this smooth quasi-static granular flow regime. 

Our approach is to test whether a straight-forward con- 
tinuum model of these smooth flow fields, based on a 
minimum of readily testable physical assumptions, can 
be made consistent with the numerical and experimen- 
tal data available for smooth quasi-static flows. In ad- 
dition to mass, linear- and angular-momentum conser- 
vation, we need to find additional relations between the 
six components of the stress tensor tJy and the state of 
the system characterized by such quantities as strain his- 
tory, packing fraction etc. Such constitutive equations 
are particularly simple for Newtonian fluids (leading to 
the Navier-Stokes equation) and elastic solids (leading 
to the equations of linear elasticity). The yielding be- 
havior of granular media illustrates that the constitutive 
equations must here take on a more complicated form. 
Granular media are athermal and dissipative — hence, 
when no external energy is supplied, grains jam into a 
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rigid, solid-like state which can sustain a finite load be- 
fore yielding [21]. Grains are made to flow by supply- 
ing an external (shear) stress to overcome this yielding 
threshold. As a result, for very slow and dense granular 
flows, the shear stresses are finite and do not approach 
zero. This complicates matters considerably. 

Our approach has two important ingredients, the de- 
tails of which can be found in Section II. Firstly, we 
are guided by the well known fact that dense grain flows 
exhibit rate independence [11]. For the velocity fields 
this means that, to good approximation, the entire ve- 
locity profile scales identically with the external driving: 
When the driving speed is doubled, the whole velocity 
field doubles. The stresses are also approximately rate 
independent, meaning that when the speed is doubled, 
the stresses stay the same. This makes the relation be- 
tween the stresses and the flow field rather special, and 
even if we could determine the full stress field, we could 
never hope to get the full velocity profile. The approach 
we take relates the stresses to certain aspects of the ge- 
ometry of the flow. This results in statements regarding 
material sheets in the flow, within which the particles 
on average only perform a collective rigid body motion 
with respect to each other. A trivial instance of such 
sheets are the layers of constant velocity present in the 
linear setups in Fig. 1 a, and b, and illustrated in Fig. 3. 
In order to say something about the actual velocities of 
such planes one would have to appeal to a sub-dominant 
dependence on shear rate. 

Secondly, when the grains are flowing, they experience 
large fluctuations [22]. Hence, we assume that if in a 
certain plane the strain rate is zero, then there will be 
no residual shear stress in this plane — if there was a 
shear stress, there would be a shear flow. Hence, all 
shear stresses are dynamically sustained, and there are 
no elastic shear stresses. Thus, we will not attempt to 
model a mixture of solid and flowing behavior as done 
in [23]. This implies that the principal- strain and stress 
directions are the same (see Section II C). 

In Section III we apply this framework to the four 
geometries depicted if Fig. 1. In the linear geometries 
(Fig. 1 a, and b), symmetry considerations directly give 
that the principal directions are constant throughout the 
system, and thus the equations are automatically closed. 
This gives the standard Mohr- Coulomb relation t = fiP, 
with the effective-friction coefficient necessarily constant 
throughout the sample. For the less symmetric geome- 
tries (Fig. 1, c and d) the local orientation of the above 
material sheets will vary throughout the cell. This allows 
us to separate the effect of constitutive assumptions re- 
garding the rates in the system and the geometry of the 
flow (see Section II B for further details) . If we maintain 
that the effective friction coefficient is constant through- 
out the sample, we find that the shear zones have in- 
finitesimal width. The standard approach of a constant 
effective-friction coefficient between shearing planes fails 
(see Section IIIB1, and especially Fig. 6). In fact, we 
then completely recover the prediction regarding the lo- 



cation of the shear zone that was derived on the basis of 
torque minimization by Unger et at in [24]. 

To capture the experimentally observed widening of 
the shear zone in the bulk, the effective-friction coeffi- 
cient has to vary throughout the shear zone. We ar- 
gue that this can only be done through a dependence 
of the effective-friction coefficient on the orientation of 
the shearing surface with respect to any bulk force (here 
gravity) . The possible origin of such an angle dependence 
is discussed in Section IV. 




FIG. 1: Geometries in which smooth, quasi-static grain flows 
occur. The velocity field is sketched with arrows, (a) Linear 
shear in absence of gravity, (b) Inclined plane flow close to 
the critical inclination angle, (c) Taylor-Couette flow with 
split bottom, (d) Linear shear over a split bottom. 



II. QUASI-STATIC GRANULAR FLOWS 

At the heart of the present development lies the pecu- 
liar fact that in order for a granular matter to support a 
shearing state, no matter how slow, a finite shear stress is 
needed. This is reminiscent of solid friction, but in stark 
contrast with the situation in Newtonian fluids. This fea- 
ture is clearly visible in the experimental results reported 
in [11, 25-27]. 

A. Explicit rate independence 

The strain-rate tensor D — (Vw -I- (Vu)')/2 plays a cen- 
tral role in ordinary fluid mechanics [28] . The use of only 
the symmetric part of the deformation-rate tensor en- 
sures that no stresses are induced by pure local rotations 
of the material (principle of material objectivity). In the 
theory of simple fluids one assumes that the knowledge 
of the complete history of D, for any material point, will 
give the stresses at that point. In the case that there are 
no memory effects this means that the stress tensor can 
be expressed as an isotropic tensor function of the strain 
rate tensor. For such functions, the first representation 
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theorem (Rivlin-Ericksen theorem [29]) states that the 
most general constitutive equation can be written as 

a = a^I + a\D + a 2 D (1) 

with a.i = ai(lD,IlD,HlD), and the invariants 

I D = tr£, I I D = trJJ 2 , III D =det£. 

As mentioned above, the granular flows that we want to 
describe are such that we have finite shear stresses even 
as the shear rate approaches zero. Thus, we can split the 
stress tensor, a, into a rate-independent part, go, and a 
rate-dependent part, gi, 

g = go+gi, (2) 

in such a way that go is not proportional to the identity 
operator (i.e. it contains shear stresses, and is hence 
not simply a hydrostatic pressure), and o\ vanishes as 
the strain rates approaches zero. The condition on go 
directly tells us that in the zero shear rate limit ol\ and/or 
a 2 must be singular. 

Theoretically there exist a flow regime in which g 
alone sets certain properties of the flow. We will refer 
to this regime as quasi-static (the precise experimental 
definition is given in Section HE). The information that 
can be extracted from the rate-independent part of the 
stress tensor will in general be of the type specifying, 
e.g., constant velocity surfaces. Due to rate indepen- 
dence, questions regarding the magnitude of the velocity 
field can not be answered by considering this limit alone, 
and neither can questions regarding the stability of any 
wide shear zones. 

We further assumes that there are only local interac- 
tions in the bulk (principle of local action) [37]. The 
implicit assumption in our approach is that our contin- 
uum description is valid, upon coarse graining over some 
small but finite length scale. 

As for the history dependence, the systems we con- 
sider are such that the flow direction is also a symme- 
try direction. For such system in steady shearing states, 
any material point will always have the same surround- 
ing flow field. Thus D does not change (up to a rotation) 
along the evolution paths of the material elements, and 
memory effects are washed out. 

The strain-rate tensor is symmetric, and is hence com- 
pletely specified by six parameters. We can choose these 
parameters as, e.g., the principal strain rates, and the 
orientation of the principal directions (specified by three 
angles). Using this parameterizations of the strain-rate 
tensor enables us to isolate the rate dependence from the 
orientational dependence. We denote the principal strain 
rates with ji, and the three angles defining the principal 
directions by 0j. The angles are to be taken with respect 
to some suitably chosen local reference direction (e.g. the 
gravitational field) . Then, the general form of the stress 
tensor is 



where the dots indicate a possible dependence on param- 
eters not directly related to the shear. Rate indepen- 
dence of the stress tensor implies invariance under the 
re-scaling v — > by, and consequently invariance under 
D — > bD 7, — ► bji). Therefore the stresses can only 
depend on the ratios of the principal-strain rates, and 
not the strain rates themselves. Hence, 

gQ = gO (7l/73, 72/73, #1,02, 03, •• •)• ( 4 ) 

To proceed with a general theory one would need to in- 
clude the full dependence on the principal-strain-rate ra- 
tios. The flows we will consider are of a limited type, 
which enables us to study the influence of the angles 
Oi without specifying the dependence on the principal- 
strain-rate ratios. We now proceed by clarifying his 
point. 

B. Shear-free sheets 

The systems we wish to consider are all such that, on 
the scale of the coarse graining, one can think of them 
as consisting of material sheets, with no internal shear, 
shearing past each other. We will here make this more 
precise and derive some important consequences. In Sec- 
tion II E we argue that the density in any flowing region is 
essentially constant, and thus mass conservation ensures 
that the flow is divergence free. This will be assumed 
already here. 

Consider a system for which it is possible to find a 
reference frame such that the velocity field is time inde- 
pendent (e.g. in the center of mass frame). We define a 
flow sheet as a surface in the flow, such that if a material 
point starts out on the surface, it stays on the surface 
throughout the time evolution of the system. If there are 
no strains within the sheet, we will refer to it as shear- 
free sheet (SFS). That is to say that the restriction of 
the strain-rate tensor to the sheet vanishes. The flows 
treated later are such that the whole shearing region can 
be divided into a collection of SFS (the SFS form a fo- 
liation of space occupied by the shear band). In any 
orthogonal and normalized basis field, with the two first 
basis vectors tangential to the SFS, the component form 
of the strain tensor is 

/ di\ 
(D) = d 2 . (5) 
V di d 2 / 

Here we have used the fact that the total flow is assumed 
to be divergence free, V • v = trD = 0. Hence, the 

principal-strain rates are 71 = and 72,3 = ±y 1 d\ + d\. 
The major advantage of considering these flows is now 
obvious: The ratio between the principal strain rates re- 
main constant throughout the system (even though d\ 
and d 2 are free to vary). Thus we can drop this depen- 
dence in stress tensor, giving 



go = go(7i,72,73,0i,02,03, ■ • ■), ( 3 ) 



go = go(0i, ■ • ■)• 



(6) 
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FIG. 2: Pictorial view of the partitioning of space into a set 
of SFS. Also illustrated is a specific instance of the SFS basis, 
and the corresponding principal-strain basis. 



The strain-rate ratios will be dropped from now on, and 
the ability to do so is crucial for the rest of the devel- 
opment. This enables us to probe the angle dependence 
alone. By a simple rotation in the SFS, the component 
form of the strain-rate tensor can be recast as 



To justify this we argue that force fluctuations are rapid 
in shearing flows. At any instance, two neighboring fluid 
elements, positioned relative to each other along any of 
the principal-strain directions, perform only a collective 
rigid-body, and a relative stretching movement. Since 
these material points are not shearing, no shear forces 
should be generated between them. If any such forces 
are present as the material points enter this no-shear 
configuration, we assume them to relax fast enough to 
be ignored. The assumption that the principal directions 
of strain and stress are aligned is also central to the flow 
rules for many of the existing continuum models of gran- 
ular flow (see [1] and references therein) . In the SFS basis 
the stress tensor takes the form 



Gro) 



SFS = 





(10) 



which again makes the connection to solid friction be- 
tween the SFS. So, the introduction of the SFS enables 
us to construct a physically relevant analogy with solid 
friction, which, as will be shown below, yields testable 
predictions. This is crucial for what remains. 
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(7) 



We will refer to the basis that realizes this compo- 
nent form of the strain-rate tensor as the SFS basis, 
{ei,e 2 ,e 3 }. Viscometric flows [28] have this form of the 
strain rate tensor, but since we will put much empha- 
sis on the physical picture offered by the SFS, we will 
continue to refer to these flows as SFS flows. This sim- 
ple form tells us that the shear between planes is always 
directed along §2- Hence, for these flows we have a pic- 
ture consistent with the SFS sliding past each other (see 
Fig. 2). For later reference the principal-strain basis (the 
basis spanned by the eigenvectors of D) {p\,P2,Pi} is 
easily seen to be given as 



Pl=§l, P2,3 = {§2 T§3)/V2, 



(8) 



in terms of the SFS basis (see Fig. 2). We now proceed 
to argue for a specific form of the stress tensor in these 
two bases. 



C. Stress relaxation 

We claim that in the principal-strain basis the stress 
tensor takes the form 
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P 1 
P 2 
OOP 3 



P i = P i (6 j ,.. 



(9) 



D. The continuity equation 

Mass conservation, and the fact that we assume the 
packing fraction to be constant throughout the shear- 
ing region, implies that the velocity field is divergence 
free. The linear-momentum continuity equation, in con- 
junction with angular-momentum continuity (and the re- 
quirement that there is no torque body couple), ensures 
that the stress tensor is symmetric, o — g} ■ This is some- 
thing we have already implicitly assumed above. The 
linear momentum continuity equation reads 



d(piQ 
dt 



v n 



(ii) 



were F is the body force, and the momentum-flux ten- 
sor is defined as II = pyv + a. As we are interested in 
quasi-static flows, we will neglect the C(|w| 2 ) term in the 
definition of the momentum-flux tensor. We will further 
only be interested in steady flows, and under these condi- 
tions the continuity equation for linear momentum takes 
the form of a force balance equation 



V-ct = F. 



(12) 



The number of additional equations needed to close such 
a system is dependent on the symmetries present, and 
will be addressed for the four geometries considered be- 
low. We now turn to dimensional analysis to determine 
what the relevant dimcnsionlcss parameters are. 
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E. Dimensional analysis and additional 
assumptions 

Obvious local parameters for the flow are the volume 
fraction <f>, the material density p mi the different local 
stresses in the SFS basis, the particle diameter a, and the 
shear rate 7. There is a further possibility that the local 
bulk-force influences the shear stress differently depend- 
ing on how it is oriented with respect to the SFS. If this is 
the case we must retain a dependence of the stress tensor 
on the orientation of the principal-strain basis with re- 
spect to the bulk force. This is encoded in the angles 0j, 
and since there is a one-to-one correspondence between 
the principal-strain basis and the SFS basis we can take 
these angle to be defined with respect to the latter in- 
stead of the former. We can now form the dimcnsionlcss 
quantities 

<P, h = t/P, v = P'/P, 6i,6 2 ,0 s , (13) 

as the shear rate is taken to zero. Since the grains are 
hard we assume the packing fraction to be independent of 
the pressure ratios, as well as the angles. Hence, we take 
the packing fraction to be constant through the quasi- 
static regime, and we drop it from the development. Ex- 
perimental and numerical justifications for this are re- 
ferred to in Section IV. The only dimensionless quantity 
that can be constructed with the strain rate is 

I=~£=- (14) 

It was shown in [11] that I is the essential parameter 
determining how the material flows. Quasi-static flow is 
to be expected for I of order 10~ 3 or less. 

As an aside we mention that for a general / in a SFS 
system (we assume 7 > for simplicity) we can write 

gi = g - go = AP'| + Ar(n/i) + &P(£>/j) 2 , (15) 

where the second equality defines the coefficients AP', 
AP, and At through (1). It is also clear that the co- 
efficients all must vanish for vanishing shear rates. For 
a SFS system only one of the fundamental invariants is 
non-zero, 

I D =0, H D =i 2 , III D =0, (16) 
and we have 




AP'( 7 ) \ 

AP( 7 ) Ar( 7 ) . (17) 
Ar(7) AP( 7 ) / 



Hence we conclude that the general form of the stress ten- 
sor is preserved even for finite shear rates. These predic- 
tions should all be possible to check by simulating these 
systems. The above forms should also be useful when 
considering the stability of these flows, and the velocity 
field on the SFS. Neither are investigated further in the 




FIG. 3: The plane shear cell, and the inclined plane geometry, 
with the SFS as well as the SFS bases indicated. 

present paper. Instead we continue to focus on the SFS 
in the quasi-static regime. 

Returning to the rate independent case, and in anal- 
ogy with solid friction, we make the additional assump- 
tion that the in-plane pressures of SFS do not affect the 
friction between the SFS. Thus we assume /i and v to 
be independent. That is, the equation of state, relating 
all the dimensionless parameters, splits into two separate 
equations 

» = n(6i), u = y{6i), (18) 

where the actual forms depend on the material. We 
therefore have four independent quantities, say P, 6\, 
62, and 63, over a three dimensional space. 

Alternatively, if we had the full velocity field of some 
suitable system, say through numerical simulations, we 
could calculate the principal-strain basis and check that 
the stress tensor has the appropriate form (9). If this 
turns out to be true, we can gain information about the 
material functions \x and v by comparing the stresses. 



III. FLOWS WITH WIDE SHEAR ZONES 

The four systems we consider (see Fig. 1) all display 
wide shear-zones and are also easily identified as SFS 
flows. We start by considering two systems with a high 
degree of symmetry, and then move on to the rather non- 
trivial split-bottom Couette geometries. 



A. Planar shear, and inclined plane geometries 

The first two geometries, to which we apply the above, 
are those of the linear shear cell without a gravitational 
field, and the inclined plane in a gravitational field (see 
Fig. 3). Due to the symmetry present in both geometries, 
we can directly identify the SFS as being parallel to the 
boundaries. Hence, e 3 is always perpendicular to these. 
The only shear present between the planes is in the di- 
rection of the velocity, so §2 points in the flow direction, 
and §1 = §2 A §3. The strain-rate tensor has the expected 
form (7) with 7 = 4^. In the considered geometries the 
principal-strain directions are constant throughout the 
sample, and thus the angles 6i are also constant. 
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The plane-shear geometry is trivial in that all the el- 
ements of the SFS basis constitute a bulk-symmetry di- 
rection. Hence, all relevant parameters must be constant 
throughout the bulk. Symmetry alone has already fixed 
the SFS planes, and appealing to subdominant depen- 
dence on /, we conclude that the shear rate is constant. 
This gives a linear velocity profile under the assumption 
that the boundaries do not break the symmetry by in- 
ducing localization of the shear zone. 

In the inclined-plane geometry the e 3 direction is no 
longer a bulk symmetry direction, and thus / is not con- 
stant along §3 . Hence we have a more complicated veloc- 
ity profile. From equation (12), or by simple force balance 
arguments between the SFS, we have \i = tan 6* (where 
9 is introduced in Fig. 3). Since the requirement on the 
effective-friction coefficient to be constant is geometri- 
cal in origin, it holds true to all orders in /. Including 
a sub-dominant dependence on / in the effective-friction 
coefficient, \i = /j,(6, 1) = tan#, thus tells us that / is con- 
stant throughout the sample. This gives the well known 
Bagnold profile [11, 17, 30]. 

Further, in this system the numerical results of [17] 
(the small inclination setups in three dimensions) show a 
linear relation between the pressures. It is also seen that 
the pressures 022 and cr 33 are very close to equal, while 
an differs substantially from the others. In the above 
treatment we have a 22 = 033 = P, in agreement with the 
numerical findings. 

In both of the above cases we have argued that the 
actual velocity profiles are set by the sub-dominant de- 
pendence on the shear rate. If true, we would expect 
strong fluctuations of the velocity around the average 
profile, something observed in both systems described 
above [11, 31]. 



is crucial to understand such slow granular flows. 



1. Rotational symmetry along the shearing direction 

The system depicted in Fig. 1 c consists of a cylindrical 
container filled with a granular material, and with a split 
bottom plate. The inner part of the container is rotated 
at an angular velocity Q, long enough for a steady state 
to be reached. The key experimental finding regards the 
spread of a wide shear band from the bottom slit up 
through the bulk to the surface. Naturally most data 
was collected for the velocity profiles at the top surface, 
as a function of the total height of the sample H . 

It was found in [18, 19] that the center position of the 
shear zone, R c , and its width, W, satisfy simple scaling 
relations as function of the layer height, H, the radial 
position of the bottom slit, R s , and the grain diameter, 
a. To good accuracy 



l-R c /R s = (H/R s ) 5 / 2 , 
Woe H^a 1 / 3 . 



(19) 



Though the experiments naturally focused on the ve- 
locity profiles at the top surface, there is also evidence 
that inside the bulk, away from the surface, the width 
of the shear zone scales with the height above the bot- 
tom z as W(z) oc z a , with a somewhere between 0.2 and 
0.4. [19, 20, 32]. 

In the natural cylindrical coordinate system, with the 
normalized basis {e r , e v , e z }, we have v = ve^, and thus 



U?)cyl 



d r oj 
d r u 8 z uj 
d?w 



ui = v/r. (20) 



B. Modified Couette geometry 

In both of the linear geometries considered above, pre- 
dicting the shape of shear zones in terms of SFS is trivial 
since symmetry guarantees that the SFS are parallel with 
the boundaries. Balancing the stresses is hence also triv- 
ial due to the special form of the stress tensor in the SFS 
basis. We now tackle the modified Couette geometry. 
Compared with the examples considered so far, this sys- 
tem has lost the symmetry in the §\ direction (along the 
SFS perpendicular to the shear; see Fig. 1). The remain- 
ing symmetry in the §2 direction is either rotational, as 
in the case of the modified Couette system (Fig. 1 c), or 
translational as in the linear system (Fig. 1 d). Though 
the loss of symmetry makes the treatment much more 
involved, it will lead to the conclusion that a constant 
effective-friction coefficient is not consistent with slow 
granular flows in general — we will find that the ap- 
propriate shape of SFS that describes the expected wide 
shear zones do not occur when we have a constant effec- 
tive friction coefficient. We suggest that a dependence of 
the friction on the local angles, as indicated in Eq. (18), 



Due to the symmetry of the problem, the surfaces of con- 
stant angular velocity are identified as the SFS. By choos- 
ing the SFS basis 



§1 = 



§2 



§3 = 



|Vw| 



|Vw| 



{d z u>e r — d r uje z ), 



(d r uie r + d z ue z ), 



(21) 



we arrive at the right form of the strain-rate tensor (7), 
with 7 = r\Vu)\/2. Due to the complicated geometry we 
need to work with the full momentum-continuity equa- 
tions (12) in order to proceed. We denote the derivative 
along the e^th direction as := e, • V. In Fig. 4 we 
have sketched a local cuboidal element of material con- 
tained between two SFS and illustrate the forces a stress 
tensor of the form Eq. (10) would give rise to. 

We introduce the angle 9 as the angle that §\ makes 
with the z-axis, K\ = gp- as the curvature of the integral 
curves of §\ (constant-w curves), and K3 = 4^ as the 
curvature of the integral curves of e 3 . Taking care of the 
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FIG. 4: Forces acting on a small element of material sand- 
wiched between two shear-free sheets (indicated). 



fact that the SFS basis forms the normalized basis vectors 
of a curved coordinate system, and using the machinery 
of tensor algebra, we can write Eq. (12) in the SFS basis 
as 



dP 



_ + (P-P')(K3-sin( 



— 3 + (ki — 2 cos#/r) r 



= —pgcos9, 

= 0, (22) 

= — pg sin 6. 



In the modified Couette setup depicted in Fig. lc, we 
only need to specify 9 to fix the SFS basis. Hence the 
relations between stress components and angles become 
of the form 



r = p(9)P, P' = v{6)P. 



(23) 



The full equations (22) coupled to Eq. (23) are too 
complicated for a full analytical treatment. We therefore 
will start from the simple assumption that the normal 
stress ratio v is equal to one (i.e., P' — P), and that p, 
is constant, i.e., independent of 9. Then we are left with 
the equations 



dP 

dx 1 
dP 

— + (m - 2 cos 9/r)P 

dP 

dx 3 



— j- = -pg cos 6, 



= 0, 



— = -pg sin6. 



(24) 



The curvature K3 has dropped out, and the first and last 
equations can be integrated to give a hydrostatic pressure 
profile P — pg(H — z). Upon substituting this into the 
second equation we conclude that the curvature of the 
constant-w curves satisfies, 



2 cos 9 sin 9 



H - z 



(25) 



To connect this formalism to the actual shapes of the 
SFS, let r(z) be the curves of constant u. Using that 



g£ = tan 6 it follows that 



sin I 

Ki -- 



_ r'(z) 



COS! 



(26) 



(l + (r'(z))2)3/2 • 

Hence using (25) and (26) we see that the curves of con- 
stant angular velocity must satisfy 

r'(z) 



r"{z) = {l + {r'(z)f) + 



H 



(27) 



This is a second order differential equation which can be 
solved numerically when supplemented with two bound- 
ary conditions. This turns out to be exactly the same 
differential equation one arrives at through minimizing 
the functional 



H 



f[r(-)} 



J dz(H 



z)r 2 {z)^l + {r'{z)f 



(28) 



As was shown by Unger and coworkers [24], this func- 
tional can be seen to describe the torque needed to 
shear an ideal cohesion-less Coulomb material under hy- 
drostatic pressure which has its infinitesimal shear zone 
at r(z). Minimizing this torque, a definite prediction 
for r(z) was obtained which qualitatively captures the 
shear zone location as measured experimentally. We re- 
fer to [19, 20, 24] for further discussion. 

However, this approach cannot result in shear zones 
with a width of the form seen in experiments, W(z) cx z a , 
with .2 < a < .5. To see this, one only has to consider 
the profiles close to the bottom. We now assume that a 
general level curve has the form [38] 



r x (z) = r Q (z) + Az a + h.o.t. 



(29) 



where ro(z) is the center curve, and A some constant 
specifying the specific level curve under consideration. 
Upon substituting this into equation (27), and consider- 
ing the lowest order in z, we conclude that a equals zero 
or one. This contradicts the experimental findings. 

From this we conclude that our assumptions that v 
and p are constant are not consistent with the wide shear 
zones observed in the modified Couette geometry. 

This curved geometry is, however, too complicated to 
study the precise role of more general v and p. We will 
therefore turn our attention to the closely related lin- 
ear split bottom shear cell, which can be obtained by 
letting the slit radius diverge, and where the rotational 
symmetry in the §2 direction is replaced by a simpler 
translational symmetry. 



2. Translation symmetry along the shearing direction 

The scaling forms (19) relate the shear zones width W, 
and location R c , to the particle size a, height H , and ra- 
dius of curvature of the slit R s . Taking the limit R s — > 00 
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SFS, constant velocity surfaces 




FIG. 5: Schematic cross section of the linear modified Couette 
geometry with the expected form of the SFS indicated. 

enables us to estimate what flow profiles can be expected 
in the linear setup shown in Fig. 1 d, even though no 
experimental or numerical data is available for such a 
system at present [39]. The width is independent of R s , 
while the shift between R s and R c should vanish in this 
limit. We therefore expect qualitatively the same widen- 
ing of the shear zone as in the Couette geometry (Fig. 1 
c), with the shear zones center remaining straight above 
the linear slit (consistent with the reflection symmetry of 
such a linear geometry), as indicated in Fig. 5. 

The equations for linear momentum conservation sim- 
plify to 

dP' 

-j-t + (P - P')k 3 = -pg cos9 
dP 

— + (P-P')K! = -pg sin9. 

Let us for the moment assume that p is a constant, 
and test whether this assumption is consistent with the 
flow profiles sketched in Fig. 5. For constant friction 
coefficient, the last two equations above can be combined 
to give 

P' Kl = pg sin e. (31) 

For profiles that bend upward when going through the 
bulk (see Fig. 5), K\ and 9 are of opposite signs. Hence, 
to satisfy Eq. (31), P' has to be negative, which is im- 
possible in cohesion-less granular materials. 

Thus we have two possible scenarios: If the effective 
friction coefficient is constant, then the system can not 
support a wide shear zone, and the shear must local- 
ize. In view of that the cylindrical geometry exhibits a 




FIG. 6: Illustration of the simple force balance arguments 
which shows that a constant p is incompatible with an upward 
bending edge of the shear zone. For details, see text. 



width of the shear zone that is apparently independent 
of the position of the bottom slit, it seems more likely 
though, that have upward bending profiles also in the 
linear geometry. Hence, the effective friction coefficient 
must decrease as we move away from the center, along 
the integral lines of §3. This is a strong statement since 
it does not rely on assuming any specific form for v 1 the 
ratio between the normal pressures: Even with normal 
stress differences we can not get a qualitatively correct 
description assuming the effective-friction coefficient to 
be constant throughout the bulk. 

The same conclusion can be reached by considering 
force balance on a cuboid element of material contained 
between two SFS as illustrated in Fig. 6 (a) and (b), 
and employing the special form of the stress tensor in 
the SFS basis (10). We start from the fact that the total 
shear forces dx l dx 2 p,P acting on the left and right edge of 
the cuboid need to balance. Now if jj, is a constant this 
implies that the normal forces on left and right of the 
cuboid equal: dx 1 dx 2 P\L = dx 1 dx 2 P\ R (Fig. 6a). The 
only additional forces acting on the cuboid are gravity 
(dx 1 dx 2 dx 3 pg) and the normal forces on top and bottom 
of the cuboid dx 2 dx 3 P'\t,b (Fig. 6b). Due to the upward 
bending of the SFS, the sum of these three terms clearly 
has a substantial component towards the right — hence 
it is impossible to balance forces in the §3 direction in this 
case. For "outward" bending SFS this problems does not 
occur as illustrated in Fig. 6 (c) and (d). 

In our case, the only way to attain force balance is if 
the normal force acting on the right face of the cuboid is 
larger than the force acting from the left: dx 1 dx 2 P'\L < 
dx 1 dx 2 P'\n- Since the shear forces still have to balance, 
this is only possible when p is not a constant — in fact 
p, has to reach its maximum along vertical SFS and then 
gradually decrease as we move outward towards increas- 
ingly slanted SFS (along the §3 direction). 

The next step is thus to include a ^-dependence in the 
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effective-friction coefficient and see if this is sufficient to 
be able to obtain shear zones of finite width. For sim- 
plicity, we keep v{6) = 1. As before, two of the mo- 
mentum continuity equations can be solved and yield a 
hydrostatic pressure profile, P — pg(H ~ z). The third 
momentum continuity equation becomes 



dln//(0) 
dx 3 



d6> 
cb 1 



sin v 
H- z 



(32) 



To get an analytically tractable problem we now consider 
a region close to the central level curve r(z) = 0. Since 
odd powers of 6 can be excluded due to the 9-^—9 
symmetry, we assume that we can expand the friction 
coefficient as 



//(#) = //,, ( 1 - ~< 



0{9 4 ). 



(33) 



Sufficiently close to the central level curve 9 is small, 
and we can, to lowest order, rewrite the derivatives ^p- 
and as 9d r + d z and — d r respectively. Hence the 
momentum conservation equation (32) can be rewritten, 
to lowest order, as 



.89 89 I 



h.o.t. 



(34) 



This differential equation can be solved by the method 
of characteristics, resulting in 



r,z) = - 



r/H 



{l + q){\- z/H)\n{\- z/H) 



+ h.o.t.. (35) 



Close to the bottom where z/H « 1, we expect the 
constant-w lines to satisfy r(z) oc z a , with an exponent 
a somewhere between 0.2 and 0.5 [19]. In this limit, (35) 
can be integrated using the fact that for any level curve 
we have r'(z) = tan 9(r(z), z). This results in 



r(z) oc z 1/(1+<?) , z/H,6<l. 



(36) 



The only sensible profiles are achieved for q > 0. This 
is in agreement with the arguments sketched in Fig. 6, 
indicating that upward bending constant-u; lines are pos- 
sible only if \x decreases with increasing \9\. Hence, the 
highest effective-friction coefficient is achieved when the 
direction of gravity lies in the tangent plane of the shear- 
ing surfaces [40]. 



IV. DISCUSSION AND CONCLUSION 

To address the question of whether continuum models 
can be made consistent with experimental and numerical 
observations of wide shear zones in slow granular flows, 
we have made a number of assumptions, which we will 
briefly recapitulate here. Trivially, we assume that the 
flow profile is smooth on a coarse grained scale. Further- 
more we assume local action (see Section II) and material 




FIG. 7: Plots of the constant-o; lines as calculated numerically 
in the small 9 approximations (33) for q = 1.5, giving an 
exponent a — 1/(1 + q) = 0.4 close to the bottom. The 
inserted graph depicts the same curves on a log-log scale, and 
with the 0.4 power indicated by the dashed line. 



objectivity (see Section II), and focus on steady states for 
which the shearing direction is a symmetry direction of 
the system — this washes out memory effects (see Sec- 
tion II). The considered flows are sufficiently symmetric, 
and time- independent, so that they can be described by 
time-independent Shear-Free Sheets. All these assump- 
tions appear rather inconspicuous. 

But there are a number of less obvious assumptions 
which deserve more attention, and for which a numerical 
test would be extremely useful. The first of these is the 
absence of elastic shear stresses in the flowing zone, due 
to rapid (on a macroscopic time scale) relaxation of force 
fluctuations (see Section II C). Clearly, in a large system, 
far away from the shear zone, elastic stresses should play 
a role, but here we only consider the actual flowing re- 
gion. It is an open question when and where such elastic 
stresses start to play a role. Secondly, we assume that 
the packing fraction is constant throughout the flowing 
region. Recent MRI measurements of the packing density 
suggest an approximately constant dilated region in the 
flowing zone [33]. Nevertheless, far away from the shear 
zone the density is observed to be different from this re- 
gion. Finally we have excluded a possible dependence of 
the effective-friction coefficient on the pressure ratio, u, 
within the SFS (see Section HE). This assumption lacks 
a strong physical argument but is made to keep the prob- 
lem tractable, and should be an important issue to check 
numerically 

Using the assumptions recapitulated above, our 
method is based on separating out those parameters of 
the strain-rate tensor that are explicitly rate dependent. 
This enables us to build a explicitly rate- independent the- 
ory, and we have shown that it is able to predict some 
of the features of the stresses seen in numerical simula- 
tions of the inclined plane geometry, as well as capturing 
the widening of the shear zone in the modified Couette 



10 



geometry. 

Through the introduction of shear-free sheets we have 
also clarified when a direct interpretation along the lines 
of solid friction can be made, and further indicated how 
far such an analogy can be stretched. Due to that the 
flow could be considered as consisting of SFS, no spe- 
cial assumptions had to be made regarding the effect of 
a variation in the principal-strain rate ratios through- 
out the sample. It was further shown that in order to 
account for the expected shape of the shear zones, the 
proportionality constants between the different pressures 
(e.g. the effective-friction coefficient) must retain a de- 
pendence on the local orientation of the flow (i.e. the 
orientation of the principal-strain basis) relative to the 
local body force — the only other probable alternative is 
that the shear zone is not wide. We speculate that the 
origin of such angle dependence is due to the competition 
between the organizational tendencies of the flow and the 
gravitational pull. The flow tends to increase the number 
of grain contacts in comprcssional directions, while de- 
crease the number in expanding directions. At the same 
time, the gravitational pull leads to an increased number 
of vertical, opposed to horizontal connections (rattlers 
always fall down). Unfortunately however, not enough is 
known about such angle dependence of the contact net- 
work in order to confirm our speculations. We suggest 
that this angle dependence as an important issue for fu- 
ture research. 

Due to the explicit rate independence of the approach, 
it can not give the complete velocity profile. In order 
to determine the complete profile one needs to include 
the sub-dominant rate dependence in the stress tensor. 



This is straight forward for some simple geometries and 
should be possible in general. Unfortunately it turns out 
to be nontrivial even for the relatively simple modified 
Couette geometry. 

Nevertheless, the intriguing fact that the experimen- 
tal shear profiles in this geometry fitted an error func- 
tion so well, provides an important benchmark for un- 
derstanding quasi-static flow. As we have discussed in 
Section IIIB2, a linear version of this experiment may 
provide important additional information. 

The present approach poses a set of well defined ques- 
tions regarding the packing fraction in the shear zone, the 
linear relationship between pressures, the simple form of 
the stress tensor in the SFS basis, and the dependence 
of the proportionality stress ratios /i and v, on the orien- 
tation of the shear planes with respect to gravity. These 
are simple basic issues which are open to investigation 
by numerical simulations, and possibly even by experi- 
ments. Clarifying these issues appears crucial for further 
development of a theory along these lines. 
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